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In  this  paper  the  central  distributional  results  of  classical  extreme 
value  theory  are  obtained,  under  appropriate  dependence  restrictions,  for 
maxima  of  continuous  parameter  stochastic  processes.  In  particular  we  prove 
the  basic  result  (here  called  Gnedenko ' s  Theorem)  concerning  the  existence 
of  just  three  types  of  non-degenerate  limiting  distributions  in  such  cases, 
and  give  necessary  and  sufficient  conditions  for  each  to  apply.  The 
development  relies,  in  part,  on  the  corresponding  known  theory  for 
stationary  sequences. 

The  general  theory  given  does  not  require  finiteness  of  the  number  of 
upcrossings  of  any  level  x.  However  when  the  number  per  unit  time  is  a.s. 
finite  and  has  a  finite  mean  y(x) ,  it  is  found  that  the  classical  criteria 
for  domains  of  attraction  apply  when  y(x)  is  used  in  lieu  of  the  tail  of  the 
marginaltfistribution  function.  The  theory  is  specialized  to  this  case  and 
applied  to  give  the  general  known  results  for  stationary  normal  processes 
(for  which  u(x)  may  or  may  not  be  finite). 

A  general  Poisson  convergence  theorem  is  given  for  high  level 
upcrossings,  together  with  its  implications  for  the  asymptotic  distributions 
of  rth  largest  local  maxima. 


2 


ley  jfagda  and  Phrases;  extreme  value  theory,  extreme  values,  stationary 

processes,  stochastic  extremes 


This  work  was  supported  by  the  Office  of  Naval  Research  under  Contract 
N00014-7S-C-0S09. 


3 

1.  Introduction. 

In  this  paper  we  shall  be  concerned  primarily  with  asymptotic 
distributional  properties  of  the  maximum 

M(T)  »  sup{?(t):  OstsT} 

of  a  continuous  parameter  stationary  process  (€(t):  t  2  0}.  A  great  deal  is 

known  about  such  properties  in  the  important  special  case  when  the  process 

is  normal  (cf.  [2],  [16]).  Our  purpose  here  is  to  delineate  the  types  of 

limiting  behavior  which  are  possible  when  the  process  is  not  necessarily 

normal,  obtaining,  in  particular,  versions  of  the  central  results  of 

classical  extreme  value  theory  which  apply  in  this  context. 

The  classical  theory  is  concerned  with  properties  of  the  maximum 

Mn  »  max(ClfC2. •••€„)  n  i.*.d.  random  variables  as  n  becomes  large. 

Central  to  the  theory  is  the  result  which  asserts  that  if  Mn  has  a 

non-degenerate  limiting  distribution  (under  linear  normalizations),  i.e.  if 

P{a  (M  -  b  )  s  x>  -*•  G(x)  for  sequences  {a  >0},{b  },  then  G  must  be  one  of 
n'nn  n  n 

only  three  general  types: 

1 

Type  I  G(x)  ■  exp(-e"x)  -°°  <  x  <  <*> 

Type  II  .  G(x)  ■  exp(-x"a)  x  >  0  a  >  0 

Type  III  G(x)  ■  exp  -(-x)a  x  <  0 

(linear  transformations  of  the  variable  x  being  permitted).  This  result, 
which  arose  from  work  of  Frechet  [S]  and  Fisher  and  Tippett  [4],  was  later 
given  a  complete  form  by  Gnedenko  [6]  and  is  here  referred  to  as  "Gnedenko's 
Theorem." 
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Gnedenko  also  obtained  necessary  and  sufficient  conditions  for  the 
doaains  of  attraction  for  each  of  the  three  limiting  types.  These  and  other 
versions  obtained  subsequently  (cf.  [7])  concern  the  rate  of  decay  of  the 
tail  l-F(x)  of  the  distribution  F  of  each  5n  as  x  increases. 

A  further  result--trivially  proved  in  the  classical  case--is  that  for 

any  sequence  {un},  r  >  0,  if  and  only  if  1  -  F(un) ~  t/n. 

This  is  sometimes  useful  in  calculation  of  the  constants  a  ,b  in  Gnedenko's 

n  n 

Theorem  (when  u  =  x/a_  +  bj . 

n  n  n 

In  more  recent  years  there  has  been  considerable  interest  in  extending 
these  and  other  results  of  the  classical  theory  to  apply  to  stationary 
sequences  which  exhibit  a  "decay  of  dependence"  which  is  not  too  slow.  In 
particular  the  early  work  of  Watson  [17]  concerning  convergence  of  P{Mn<un> 
applied  under  m-dependence ,  Loynes  [14]  proved  Gnedenko’s  Theorem  under 
strong  mixing  assumptions,  and  Berman  [1]  obtained  detailed  results  for 
normal  sequences  under  a  mild  condition  involving  correlation  decay.  More 
recently  we  have  obtained  a  theory  (cf.  [9])  involving  weak  "distributional 
mixing"  conditions,  which  unifies  these  results  and  provides  a  rather 
satisfying  extension  of  the  classical  distributional  theory  to  include 
stationary  sequences. 

It  is  not  too  surprising  that  such  an  extension  is  possible  for 
stationary  sequences,  at  least  under  suitable  dependence  restrictions.  What 
may  seem  surprising  at  first  sight  is  that  a  corresponding  theory  is 
possible  for  continuous  parameter  stationary  processes.  However  this 
becomes  intuitively  clear  by  recognizing  that  the  maximum  up  to  time  n,  say, 
is  just  the  maximum  of  n  random  variables--the  "submaxima"  in  the  fixed 
intervals  (i-1  ,  i),  1  s  i  i  n.  Our  procedure  will  be,  in  fact,  to  use  the 


existing  theory  for  stationary  sequences  by  aeans  of  (a  slightly  modified 
version  of)  this  precise  approach.  The  sequence  results  which  will  be 
needed  are  stated  in  Section  2. 

In  Section  3  we  will  obtain  Gnedenko's  Theorem  for  continuous  parameter 
stationary  processes,  showing  under  appropriate  conditions  that  if 

P{aT(M(T)-bT)  s  x}  G(x)  »s  T  * 

for  some  constants  a^,>0,bT,  then  G  must  be  one  of  the  extreme  value  forms. 

In  Section  4  we  obtain  a  related  result--again  extending  a  classical 
theorem--to  give  necessary  and  sufficient  conditions  for  the  convergence  of 
P{M(T)$Uj.}  for  sequences  not  necessarily  of  the  form  Uj.  ■  x/aT  ♦  bT 
implicit  in  Gnedenko's  Theorem. 

As  a  corollary  of  this  result  we  obtain  necessary  and  sufficient 
criteria  for  the  domains  of  attraction  occurring  in  Gnedenko's  Theorem.  In 
the  classical  i.i.d.  sequence  case,  the  criteria  for  domains  of  attraction 
involve  the  rate  of  decay  of  the  marginal  distribution  l-F(x)  as  x 
increases.  For  the  present  case  the  very  same  criteria  apply,  provided 
l-F(x)  is  replaced  by  another  function  ip(x).  For  processes  whose  mean 
number  y(x)  of  upcrossings  of  any  level  x  is  finite,  the  function  tp(x)  is 
precisely  y(x),  a  readily  calculated  quantity. 

The  general  theory  will  not  require  that  the  mean  number  of  upcrossings 
of  a  level  per  unit  time  be  finite,  and  accordingly  will  include  the  class 
of  stationary  Gaussian  processes  with  covariances  of  the  form 
r(t)  -  1  -  C|t|°  ♦  o(|t|°)  as  T  0  for  0  <  a  <  2.  In  Section  5  we  consider 
such  processes,  as  well  as  (possibly  non-Gaussian)  cases  for  which  the  mean 
number  of  upcrossings  per  unit  time  ia  finite.  Finally  in  Section  6  we  note 


the  general  Poisson  limit  for  the  point  processes  of  uperossings  of 
increasingly  high  levels  and  its  implications  regarding  limit  theorems  for 
the  distribution  of  the  r**1  largest  local  maximum  of  £(t)  in  0  s  t  s  T. 


2.  Two  results  for  stationary  sequences. 

As  noted,  our  development  of  extremal  theory  for  stationary  processes 
will  rely  in  part  on  the  existing  sequence  theory.  Specifically  we  shall 
require  the  following  definitions  and  results  (which  may  be  found  e.g.  in 


Let  {£_ }  be  a  stationary  sequence  and  write  F.  .  (x,...x  )  for  the 

at  la  I  »  •  X  _  X  ll 

1  n 

joint  distribution  function  of  £.  .  For  brevity  write  also  F.  .  (u) 

*i  Xq  •  •  •  x 

l  xi  . ,  —  1  n 

to  denote  F.  .  (u,u...u)  *  P{£,  Su  Su>.  If  {u  }  is  a  sequence  of 

1l-,'1n  X1  ln  n 

real  constants,  we  say  that  the  sequence  {£n>  eaHsfiee  the  (dependence) 

condition  D(un)  if  for  each  n,  1  s  ij<  ij. .  .<  ip<  jj. . . <  jp,  £  n  ,  jj.  i  i  l, 
(2.1)  |F  (u  )  -  F  (u  )F.  .  (u  )  |  s  a 

‘j*  •  pJ  i*  •  •  J  pi  "  M’”  p  n  ^’“‘■'p'  "  "»*' 


where 


(2.2) 


a  .  ■*  0  for  some  sequence  l  »  0  as  n 
n,tn  n 


Note  that  a  g  can  (and  will)  be  taken  to  be  decreasing  in  l  for  each  n 
by  simply  replacing  it  by  the  smallest  value  it  can  take  to  make  (2.1)  hold 
(i.e.  the  maximum  value  of  the  left-hand  side  of  (2.1)  over  all  allowable 
sets  of  integers  i^...^,  jj...jp,  .  Note  also  that  (2.2)  may  then  be  shown 
equivalent  to  the  condition  (cf.  [12]  for  proof) 


(2.3) 


a„  +  0  as  n  •  for  each  X  >  0  . 
n,nA 
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The  condition  D(un)  indicates  a  degree  of  "approximate  independence"  of 
members  of  the  sequence  separated  by  increasing  distances.  However  this 
condition,  which  we  refer  to  as  "distributional  mixing,"  is  clearly 
potentially  far  less  restrictive  than,  for  example,  "strong  mixing."  In  the 
case  of  normal  sequences,  it  is  in  fact  satisfied  when  the  covariance 
sequence  {rR}  tends  to  zero  even  just  fast  enough  so  that  rn  logn  0. 

The  following  result  is  basic  to  the  sequence  theory  and  will  be 
required  in  later  sections. 

Lemma  2.1.  Let  {^}  be  a  stationary  sequence  satisfying  D{un>  for  a  given 
sequence  {ur}  of  constants  and  write  Mn  ■  maxfCj  ,S2*  •  Then  for  any 

integer  k  2  1  (writing  []  to  denote  integer  part), 

P^MnSun*  "  P^M  [n/k]Sun^  *  0  08  n  *  *  • 

This  lemma  indicates  a  degree  of  independence  between  the  [n/k]  maxima 
when  the  first  n  integers  are  divided  into  k  groups.  We  shall  also  need  the 
sequence  form  of  Gnedenko's  Theorem,  which  is  given  (e.g.  in  [10])  as  follows: 

Theorem  2.2.  Let  {£n>  be  a  stationary  sequence  such  that  Mn  *  maxfCj,^.  •  -Cn) 

satisfies  ?{•„(**„" bn)  s  x)  -*•  G(x)  «  n  +  «  for  some  non-degenerate  d.f.  G 

and  constants  {a  >  0},{b  }  .  Suppose  that  D(u  )  holds  for  all  u  of  the  form 
n  n  n  n 

x/a_  +  b_,  -eo  <  x  <  00 .  Then  G  is  one  of  the  three  extreme  value 
iin 

distributional  types. 

The  other  classical  result  quoted— concerning  convergence  of  P(Mnsun) 
for  arbitrary  sequences  {un>— Is  also  important  and  holds  under  appropriate 
conditions  for  stationary  sequences  {€„)•  This  will  not  be  discussed  here 
since  the  corresponding  continuous  parameter  result  will  be  independently 
derived. 
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5.  Gnedenko's  Theorem  for  stationary  processes. 

As  indicated  above,  it  will  be  convenient  to  relate  the  maxima  M(T)  of 
the  continuous  parameter  stationary  process  £(t)  to  the  maximum  of  n  terms 
of  a  sequence  of  "subaaxima."  Specifically  if  h  >  0  we  write 

(3.1)  sup{£(t):  (i-l)hstsih) 
so  that  for  n  ■  1,2,3...  , 

(3.2)  M(nh)  ■  max(C1,C2,...Sn)  . 

The  following  preliminary  fora  of  Gnedenko's  Theorem  (involving 
conditions  on  the  C-sequence)  is  immediate. 

Theorem  3.1.  Suppose  that  fop  some  families  of  constants  {aT>0),{bT>  we 
home 

(3.3)  P{aT(M(T)-bT)  £  x}  G(x)  as  T  » 

for  seme  non-degenerate  G,  and  that  the  { )  sequence  defined  by  (3.1) 
satisfies  D(un)  whenever  un  «  x/a^  +  b^  for  some  fixed  h  >  0  and  all  real 
x.  Then  G  is  one  of  the  three  extreme  value  types. 

Proof.  Since  (3.3)  holds  in  particular  »s  T  +  »  through  values  nh  and  the 
sequence  is  clearly  stationary,  the  result  follows  by  replacing  £n  by  Cn 
in  Theorem  2.2  and  using  (3.2).  □ 

Corollary  3.2.  The  result  holds  in  particular  if  the  D(un)  conditions  are 
replaced  by  the  assumption  that  (£(t)}  is  strongly  mixing.  For  then  the 
sequence  (cn>  is  strongly  mixing  and  satisfies  D(un). 


0 
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We  now  introduce  the  continuous  analog  of  the  condition  D(un) ,  stated 

in  tens  of  the  finite  diaensional  distribution  functions  F«,  of  £(t) 

tl***tn 

(again  writing  F  (u)  for  F  (u...u). 

i  n  in 

The  oondition  D  (u^.)  will  be  said  to  hold  for  the  process  £(t)  cmd  the 
family  of  constants  (uj.:  T>0},  with  respeot  to  a  family  {q^. 0},  if  for 
any  points  Sj<  s2« .  .<  s^<  tj. .  .<  tp,  belonging  to  (kqT:  Oskq^.sT)  and 
satisfying  t ^  i  t,  we  have 


(3.4)  |F 


V 


’“pV 


t  <V  ' 

p' 


Vv 


t 

•p* 


s  a 


T,Y 


where 


-*■  0  for  some  sequence  yt 


O(T)  or,  equivalently,  where 


(3.5) 


“t.at  *  ®  “  T  *  " 


for  each  X  >  0. 

The  D(un)  condition  for  (Cn)  required  in  Theorem  3.1  will  now  be 
related  to  D^u^)  by  approximating  crossings  and  extremes  of  the  continuous 
parameter  process,  by  corresponding  quantities  for  a  sampled  version.  To 
achieve  the  approximation  we  require  two  conditions  involving  the  maximum  of 
€(t)  in  fixed  and  in  very  small  time  intervals.  These  conditions  are  given 
here  in  a  form  which  applies  very  generally--readily  verifiable  sufficient 
conditions  for  important  cases  are  given  in  Section  5. 

Specifically  we  suppose  that  there  is  a  function  t{/(u)  such  that,  for 
h  >  0, 

(3.6)  limjwp  *  1  . 


♦ 
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and  that  for  each  a  >  0,  there  is  a  family  of  constants  q  *  q&(u)  +  0  as 
u  <*>  such  that 


(3.7) 


lim  sup 

lf400 


PU(0)<u 


q'l'(u) 


-*■  0  as  a  0 


Note  that  Equation  (3.6)  specifies  an  asymptotic  upper  bound  for  the 
tail  distribution  of  the  maximum  in  a  fixed  interval,  whereas  (3.7)  limits 
the  probability  that  the  maximum  in  a  short  interval  exceeds  u,  but  the 
process  itself  is  less  than  u  at  both  endpoints.  The  following  result  now 
enables  us  to  approximate  the  maximum  in  an  interval  of  length  h  by  the 
maximum  at  discrete  points  in  that  interval. 

Lemma  3,3. 

(i)  If  (3.6)  holds,  then  P{M(q)>u)  =  oip(u)  ae  u  ■+•  »  for  any 
q  =  q(u)  -*■  0.  Also  P(5(0)>u)  =  oiJ>(u). 

(ii)  If  (3.6)  and  (3.7)  both  hold ,  and  I  te  an  interval  of  length  h, 
then  there  are  constants  X  such  that 

A 

(3.8)  0  s  lim  sup[P{5(jq)s  u  ,  jqel)  -  P{M(I)su)]/i(»(u)  s  X  -*•  0  as  a  -*■  0  , 

u^«>  a 

where  q  *  q  (u)  is  as  in  (3.7),  the  convergence  being  uniform  in  all  intervals 

A 

of  this  fixed  length  h. 

Proof.  If  (3.6)  holds  and  q  +  0  as  u  +  «»  then  for  any  fixed  h  >  0,  q  is 
eventually  smaller  than  h  and  P(M(q)>u}  s  P{M(h)>u),  so  that 

lim  sup  P{M(q)>u}/ip(u)  s  lim  sup  P(M(h)>  u}/tp(u)  s  h  by  (3.6)  , 

U-*°°  u-*« 

from  which  it  follows  that  P{M(q)>u}/iJ/(u)  0,  as  stated.  The  remaining 

statement  of  (i)  also  follows  since  P{£(0)>u}  i  P{M(q)>u). 
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Suppose  now  that  (3.6)  and  (3.7)  both  hold  and  that  I  is  an  interval  of 
fixed  length  h.  The  interval  I  consists  of  no  more  than  h/q  subintervals  of 
the  form  ((j-l)q  ,  jq),  together  with  (possibly)  a  shorter  interval  at  each 
end.  The  difference  in  probabilities  in  (3.8)  is  clearly  non-negative  and 
(using  stationarity)  dominated  by 

K  „  *  irPtt(0)<u,e(q)<ufM(q)>u}  ♦  2P{M(q)>u}  . 

a,u  q 

The  desired  result  (ii)  now  follows  from  (3.7)  and  (i)  by  writing 

X  *  lim  sup  X  .  □ 

a  a  «u 

*  ir**0  » 

It  is  now  relatively  straightforward  to  relate  D(un)  for  the  sequence 
{£n>  to  the  condition  for  the  process  £(t),  as  the  following  lemma 

shows.  (In  this  we  use  the  (potentially  ambiguous)  notation  Dfu^)  to  mean 
D(v„)  «ith  »„.  unh.) 

Lemma  3.4.  Suppose  that  (3.6)  holds  with  some  function  i|>(u)  let  (qa(u)} 
be  a  family  of  constants  for  each  a  >  0  with  qa(u)>0,  q&(u)-»-0  as  u  -*•  °°, 
and  such  that  (3.7)  holds.  If  D^u^,)  is  satisfied  with  respect  to  the 
family  qT  =  q&(uT)  for  each  a  >  0,  and  TiKUj)  is  bounded ,  then  the  sequence 
{^n>  defined  by  (3.2)  satisfies  D(unh)  for  h  >  0. 

Proof.  For  a  given  n,  let  i^<  ij. .  .<  ip<  .  .<  jp,<n,  jjj  -  ip  *  £ .  Write 

Ir“  [(ir-i)h  »  irhl  >  Js  *  C (j *  ish3’  For  brevity  let  q  denote  one  of  the 

families  {q  (•)}  and 
& 

Aq  *  n  {£(jq)su  jqclr},  AUi/Unh> 

n  r#l  r* i  r 

p»  p' 

Bn-  n  U(jq)Su  .  ,  jqeJ  )  ,  B-  n  U.  Su  .}  . 
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It  follows  in  an  obvious  way  from  Lemma  3.3  that 

0  s  lim  sup{P(A  nB  )  -  P(AnB) }  i  lim  sup(p«-p')<|>(unh)Xa 
n-*«  **  **  n-*» 

s  ”«u„lA  1 

for  some  constant  K  (since  nhij)(unh)  is  bounded)  and  where  X.,  0  as  a  •+•  0. 

Similarly 

lim  sup|P(Aq)-P(A)|iKXa  ,  lim  sup|P(Bq)-P(B)  | as  KXft  . 

Now 

|P(AnB)  -  P(A)P(B) |  S  |P(AnB)  -  P(AqnBq) |  ♦  |P(AqnBq)  -  P(Aq)P(Bq)| 

(3.9)  ♦  P(Aq)|P(Bq)  -  P(B) |  ♦  P(B)|P(Aq)  -  P(A)| 

■  *  I"<VV  •  P(VP(V  I 

where  lim  sup  Rn  a  s  3KX#. 

n-*«  ’ 

Since  the  largest  jq  in  any  IT  is  at  most  i^.h,  and  the  smallest  in  any 

J  is  at  least  (j-l)h,  their  difference  is  at  least  (*-l)h.  Also  the 
1 

largest  jq  in  Jp,  does  not  exceed  jp,h  s  nh  so  that  from  (3.4)  and  (3.9) 

(3.10)  |P(AnB)  P(A)P(»|  * 

(in  which  the  dependence  of  aT  ^  on  a  is  explicitly  indicated).  Write  now 
a*  t  ■  inffRjj  #  ♦  «nh\£-l)h^‘  Since  the  left-hand  side  of  (3.9)  does  not 

depend  on  a  we  have 


|P(AnB)  -  P(A)P(B) |  s  a*  ^  , 


IS 


which  is  precisely  the  desired  conclusion  of  the  lemma,  provided  we  can  show 
that  lim  a*  .  ■  0  for  any  X  >  0  (cf.  (2.3)).  But  for  any  a  >  0 

Mm  "» An 


<  D  4.  /,^)  <  D  *  rt(^) 

n,An  n,a  nh,(Xn-l)h  n,a  nh.JjXnh 


when  n  is  sufficiently  large  (since  decreases  in  £),  and  hence  by  (3.5) 


lim  sup  a*  ,  s  3KX_  , 
n-H»  n,*n  a 


and  since  a  is  arbitrary  and  A„  0  as  a  0,  it  follows  that  a*  .  0  as 

a  n.  An 

desired.  Q 

The  general  continuous  version  of  Gnedenko's  Theorem  is  now  readily 
restated  in  terms  of  conditions  on  £(t)  itself. 

Theorem  3.5.  With  the  above  notation  for  the  stationary  prooeas  5(t) 
satisfying  (3.6)  for  some  function  suppose  that,  for  some  families  of 
constants  {a^oMhj,}, 


P{aT(M(T)-bT)  s:  x)  -*■  G(x) 

for  a  non-degenerate  G.  Suppose  that  Tip(u^)  is  bounded  and  Dc(Uj.)  holds  for 
Oj.  ■  x/Sj  *  bT  for  each  real  x,  with  respect  to  families  of  constants 
(qa(u)}  satisfying  (3.7).  Then  G  is  one  of  the  three  extreme  value 
distributional  types. 


Proof.  This  follows  at  once  fro*  Theorem  3.1  and  Lemma  3.4. 


0 
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As  noted  the  conditions  of  this  theorem  are  of  a  general  kind,  and  more 
specific  sufficient  conditions  will  be  given  in  the  applications  in 
Section  5. 

4.  Convergence  of  PtMCnsuj,}. 

Gnedenko's  Theorem  involved  consideration  of  P{a.j,{M(T)-bT}  £  x},  which 
may  be  rewritten  as  P{M(T)£u,p}  with  ttj,  *  a^,1  x  +  bT.  We  turn  now  to  the 
question  of  convergence  of  PtMCnsu^.}  as  T  -*•  »  for  families  u^.  which  are 
not  necessarily  linear  functions  of  a  parameter  x.  (This  is  analogous  to 
the  convergence  of  p(MnSun)  for  sequences,  of  course.)  These  results  are 
of  interest  in  their  own  right,  but  also  since  they  make  it  possible  to 
simply  modify  the  classical  criteria  for  domains  of  attraction  to  the  three 
limiting  distributions,  to  apply  in  this  continuous  parameter  context. 

The  discussion  will  be  carried  out  in  terms  of  so-called  "e-upcrossings” 
of  a  level  by  the  stationary  process — a  concept  originally  introduced  by 
Pickands  [15]  to  deal  with  extremes  of  processes  whose  sample  functions  were 
so  irregular  that  the  "ordinary"  upcrossings  could  be  infinite  in  number  in 
a  finite  interval.  (Here  we  make  essential  use  of  this  concept  whether  the 
process  is  irregular  or  not.) 

Briefly,  if  e  >  0,  C(t)  is  said  to  have  an  e-upcrossing  of  u  at  a  point 
tQ  if  5(t)  s  u  for  all  t  in  the  interval  (tQ-e,  tQ),  but  £(t)  >  u  for  some 
point  t  e  (t^  ,  tQ+n)  for  each  H  >  0.  Since  the  interval  (tg-e,tg)  contains 
no  upcrossings,  the  number  of  e-upcrossings  in  a  unit  interval  does  not 
exceed  1/e  .  We  write  N  (t)  ,  N  (I)  for  the  number  of  e-upcrossings  in 

v»  p  U  W|U 

the  intervals  (0,t),I  respectively  and  y  •  EN  (1)  so  that 

t  )U  C|U 

EN  .(f)  *  ty  .  The  following  small  result  indicates  some  connections 
between  e-upcrossings  and  maxima. 
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La—  4.1. 

(i)  For  h  >  0,  P{M(h)>u}  2  hvh  u  . 

(H)  hy.  a  P(M(2h)>u)  -  P{M(h)>u). 
n,u 

(Hi)  If  (3.6)  holds  (i.e.  lim  sup  P{M(h)>u)/(htJ»(u))  £  1  for  some  ip(u), 

u-*“ 

h  >  0)  then  lim  sup  y.  /( |>(u)  s  1. 

U-K»  n,U 

(iv)  If 


(4.1)  P{M(h)>u)  ~  htji(u)  «iu  +  «  for  0  <  h  <  hQ 

then  y  ~  ^(u)  for  all  ( sufficiently  email)  e  >  0. 

£,U 


Proof.  Since  clearly  u(h)  is  either  zero  or  one,  we  have 

hyh>u  -  ENh,u(h)  *  p{Nh,u(h)“1}  *  p{M<h>>u} 

so  that  (i)  follows  at  once.  To  prove  (ii)  we  note  that 
P{M(2h)>  u)  s  P(M(h)>u}  +  P{N.  fh,2h)2s  1} 

*i|U 

•  P{M(h)>u)  +  hyh  u  . 

If  (3.6)  holds  (Hi)  follows  at  once  from  (i) . 

Finally,  if  (4.1)  holds,  the  conclusion  of  (iv)  follows  from  (Hi)  and 
the  inequality  (ii) ,  which  gives 


lim  inf  y  „/iKu)  *  lim 

u-k»  e,u  u-h» 


f 


P{M(2e)>u>  P(M(e)>u) 

c^(u)  "  eiJu) 


] 


□ 


Our  main  purpose  is  to  demonstrate  the  equivalence  of  the  relations 
P{M(h)>uT)  ~  t/T  and  P{M(T)SuT>  e"T  under  appropriate  conditions.  The 
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following  condition  will  be  referred  to  as  D'(uT),  and  is  analogous  to  D' 

C  I 

conditions  required  for  similar  purposes  for  sequences  (cf.  [10]). 

If  {tij.}  is  a  given  family  of  constants  the  condition  0*(Uf)  will  be 
said  to  hold  (for  the  process  (£(t)}  satisfying  (4.1))  if 


(4.2) 


lim  sup  T|yeT  -  iKi^)  j  ♦  0  as  c  -*■  0 


We  now  state  and  prove  the  first  part  of  the  desired  equivalence. 

Theorem  4.2.  Suppose  that  (4.1)  holds  for  some  function  ij»  and  let  {u^,}  be  a 
family  of  constants  such  that  DC(U|.)  holds  (with  respect  to  a  family  Cq) 
satisfying  (3.7)),  and  that  D*^)  holds.  Then 


implies 


TiKUj.)  T  >  0 


PlMCHSu,.}  ♦  e' 


Proof .  Let  0  <  h  <  h^  (cf.  (4.1)),  and  let  n,k  be  integers,  writing 
n'  ■  [n/k].  By  Lemma  3.4,  the  sequence  of  "submaxima"  (cn)  defined  by 
(3.1)  satisfies  D(unh)  and  hence,  from  Lemma  2.1, 


(4.S)  P{M(nh)sunh)  -  P*{M(n'h)sunh}  0  as  n  ♦  -  . 


Now  writing  u  .  -  u,  Lemma  4.1  (i)  gives 


n'hVh,u  *  p<#(n*h)>u)  «  n'P{M(h)>u)  ~  n  •!*(«) 


so  that 
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(4.6)  1  -  n'h*(u)(Uo(l))  i  P(M(n'h)su)  s  1  -  n'l»Vh,u  . 

But  n'h^(u)  ■  nhi(i(unh)/k  r/k  by  (4.3).  Further 


n'h,u 


Ttunh 

T,unh 


*(,W 


so  that  letting  n  -►  •  in  (4.6)  and  using  D^(u^.)  we  have 

1  -  t/k  s  lim  inf  P(M(n'h)su)  s  lia  sup  P{M(n'h)su) 

n-«a  n-*o> 

si-  t/k  ♦  o(l/k)  . 

By  taking  kth  powers  and  using  (4.5)  we  see  that 

(l-T/k)k  s  lia  inf  P{M(nh)Su)  s  lin  sup  P{M(nh)su}  s  (1-  t/k  +  o(l/k)) 
and  hence,  letting  k  -*•  ®,  that 


(4.7)  P{M(nh)sunh)  ♦  e  . 

Now  if  n  is  chosen  so  that  nh  s  T  <  (n+l)h,  and  if  uRh  s  Uj., 
P{M(nh)SuT)  ■  P{M(nh)sunll>  ♦  Piu^sM^^Su^.} 


where  the  last  term  does  not  exceed 


nPtunhSM(h)<uT}  B  ntp{MW>unh}  "  p^(h)>ur}] 

■  [g  (!♦•(!))  -  ^  (1+0(1))] 


by  (4.1)  and  (4.3).  Since  nh  ~  T  this  clearly  tends  to  zero.  A 
corresponding  calculation  where  >  u^.  thus  shows  from  (4.7)  that 

P(M(nh)suT>  -►  e"T 

with  T  »  [n/h].  Finally 

POKT)*^)  s  P{M(nh)suT) 

s  ♦  P(H(nh)suT<M(T)}  . 

But  the  event  M(nh)  z  uT  <  M(T)  implies  that  the  maximum  in  the  interval 
[nh, (n+l)h]  exceeds  Oj.,  which  has  probability  P{M(h)>u_},  giving 

P(M(nh)suT)  -  P(M(h)>uT)  S  P{M(T)sutJ  s  P{M(nh)SuT)  , 

from  which  (4.4)  follows  since  P{M(h)>Uj,}  ~  hx/T  -►0.  □ 

In  our  treatment  of  the  converse  result  it  will  be  convenient  to  use 
the  innocuous  further  assumption 

(4.8)  ~  ^^u[x/h]h^  as  n  * 

(for  some  given  h  >  0).  This  assumption  is  possibly  dispensable  but 

L 

certainly  commonly  holds  (e.g.  when  \J>(uT)  ~  (2  log  T)  for  stationary 
normal  processes)  and,  of  course,  always  holds  if  the  function  Uj.  is 
replaced  by  the  step  function  u[T/h]h’  constant  between  consecutive  points  nh. 

The  first  step  of  the  derivation  exhibits  approximate  independence  of 
maxima  in  disjoint  intervals. 
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La— »  4.3.  Suppooe  that  (4.1)  hols1*  for  ocmm  ^  and  lot  {i^}  bo  a  family  of 
oonotanto  (oatiofying  (4.4)  for  oomo  h  >  0)  ruoh  that  Dc(u^,)  holds  with 
roopoot  to  a  family  (q)  oatiofying  (S.7)  and  ruoh  that  TiKttj.)  is  bounded. 
Let 

(4.9)  P{M(T)iuT)  -►  e"T 
for  some  x  >  0,  Then  for  k  •  1,2... 

(4.10)  P{M(T/k)Su_}  -*•  e"T/k  ao  T  ■*■<*>  . 

Proof.  As  in  the  previous  result  the  assumptions  surrounding  Dc  (u^.)  imply 
that 

(4.11)  P{M(nh)Su}  -  Pk{M(n'h)Su}  ■*  0  as  n  -*•  ® 

where  n'  *  [n/k]  and  u  -  u  Now  if  n  »  [T/h]  it.  is  readily  checked  that 
[n/k]h  s  T/k  <  ([n/k]+l)h,  so  that 

P{M(T/k)Su}  s:  P{M(n'h)su)  -*>  e‘T/k 

by  (4.5)  (which  holds  here  by  the  same  argument  as  in  Theorem  4.2),  and 
(4.9)  with  T  ■  nh.  But  also 

P{M(T/k)Su}  2  P{M((n'+l)h)Su> 

-  P{M(n'h)su)  -  P{M(n'h)s  u<M((n,+l)h)  } 

2  P{M(n*h)su)  -  P(M(h)>u)  . 

The  first  term  on  the  right  tends  to  e'T/k  ,  and  the  second  is  asymptotically 
equivalent  to  h*(unh)  -  n_1(nh«(unh))  0  since  *np(uT)  is  bounded.  Hence 

we  have 
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(4.12)  P{M(T/k)Sunh)  •*  e‘T/k  . 

But  it  follows  simply  that  u^  may  be  replaced  by  Up  in  (4.12)  to  give 
the  desired  result  since  if,  for  example,  un^  s  Uj,  we  have 

0  s  P{M(T/k)SuT)  -  P{M(T/k)sunh) 

(4.13)  =  P*unh<M(T/k)SuT} 

s  (n/k)P{unh<  M(h) s u^j.}  *  P{M(h)>unh> 

since  if  the  maximum  in  (0  ,  T/k)  lies  between  unh  and  Uj.  this  must  also 
occur  in  one  of  the  first  n'  intervals  ((i-l)h  ,  ih)  or  in  (n'h  ,  T/k).  The 
first  term  of  (4.13)  is  readily  seen  to  be 

(nh/k)  [<|/(unh)  (1*0(1))  -  (Kit,.)  (l+o(l))]  , 

which  is  easily  seen  to  tend  to  zero  by  (4.8)  since  nh\p(unh)  is  bounded. 
Boundedness  of  nhKu^)  also  implies  that  the  second  term  of  (4.13)  tends 
to  zero.  A  corresponding  calculation  applies  for  unh  2  uT  so  that 
PtMCT/kjSitp}  -  PtMCT/kJsu^}  -*■  0,  giving  the  desired  result.  □ 

Lemma  4.4.  Under  the  seme  aeexmptione  ae  in  Lemma  4.3  we  have 

i  lim  sup|P{M(eT)suT>  -  e"eT|  *  0  ae  t  ♦  0  . 

T-*» 

Proof.  Choose  the  integer  k  depending  on  e  >  0,  so  that  s  e  <  y  . 
Then 


PtMa/^SUy}  i  P{M(eT)iuT>  s  P{M(T/(k*l))suT)  . 
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By  subtracting  e”ET  from  each  tern,  and  using  the  facts  from  Leans  4.3  that 
P{M(T/k)suT}-*e‘T/k  ,  P{M(T/(k+l))suT}+e"T/(k+1),  we  see  that 

lia  sup|P{M(eT)s  u_}  -  e_eT)  s  aax[|e'T/(k+1)-  e"eT|  ,  |e‘T/k  -  e”ET|]  . 

But 

i|,-t/(»*l)  .  ,-et|  t  i(1  .  e-r(e  -  1/Ck.l))| 

.  I  [i  .  ,-i(e-t(l.cO)))| 

which  tends  to  zero  as  e  -*•  0  .  Similarly  i|e~T^k  -  e"eT|  +  0  so  that  the 
desired  result  follows.  □ 

The  next  lemma  gives  a  conclusion  which  is  interesting  in  itself  and 
from  which  the  main  result  will  follow  immediately. 

Lemna  4.5.  Again  suppose  that  the  conditions  of  Lemma  4.3  hold.  Then 

(4.14)  lim  sup | Ty  _  -  x|  0  as  e  -*■  0  . 

T-Ka  e*»ur 

Proof.  By  Leiana  4.1  (if), 

P{M(2cT)>ut)  s  P{M(eT)>uT)  ♦  eTueT 

so  that 

TM-t  „  *  T  2  £  (PCMfeTJsu^.}  -  e'ET] 

*  T 

♦  [e“2eT  -  P{M(2eT)suT)J  ♦  I[e_ET  -  e*2ET  -  ex]  , 
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giving 

lin  inf [Ty  _  -  x]  *  i[e-£T  -  e_2eT  -  ex] 

T-KO  e,»“T  6 

-  i  lim  supIPfMCeTJSi^}  -  e"eT| 

T-mo 

-  |  lim  sup|P{M(2eTjsuT>  -  e"2€T|  . 

T-*» 

The  latter  two  terms  tend  to  zero  as  £  ■*  0  by  Lemma  4.4,  so  that 
lim  inf[Tu  _  -  x]  &  a  -*•  0  as  e  -*•  0  . 

f-»oo  *'1,U'T  e 

Similarly  from  the  inequality  £Ty| 
follows  that 

lim  sup[TucT  -  T]  s  b 

T-*«  61  »“t  e 

where  b  0  as  e  0.  Since  if  lim  inf  8  £  X  and  lim  sup  8  s  A  it  is 
t  n  *  n 

easily  shown  that  lim  sup|8n|  s  max(J A| , (A|) ,  we  have 
lin  sup|Ty  -  x|  s  max(|a  |,|b  |)  , 

C 1  *  UT  **  fc 

which  tends  to  zero  as  e  0,  giving  the  conclusion  of  the  lemma.  □ 

It  will  be  noted  that  (4.14)  is  very  similar  to  the  condition  ^(Uj) 
and  follows  as  a  conclusion  from  the  assumption  PtMCOsu^,}  e"T  under 
appropriate  conditions.  For  this  we  do  not  require  that  D^fu^.)  hold.  If 
we  now  do  assume  that  0^(1^)  holds  we  immediately  obtain  the  main  converse 


£  P(M(eT)>uT}  (Lemma  (4.1)  (i))  it 


result. 
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Theorem  4.6.  Suppose  that  (4.1)  holds  for  some  t  and  let  (u^J  be  a  family 
of  constants  (satisfying  (4.8)  for  some  h  >  0)  suoh  that  Dc (u^.)  holds  with 
respeot  to  a  family  (q)  satisfying  (3.7).  Suppose  also  that  D£(Uj.)  holds. 
Then  (4.4)  implies  (4.3). 


Proof.  Since  Ty„„,  ..  s  T/(eT)  ■  1/e  it  is  implicit  in  the  assumption 
"  »u<j« 

(tUp)  that  Tty(Uj.)  is  bounded.  Thus  the  conditions  of  the  previous  lemma 
are  satisfied  if  (4.4)  holds  and  hence 


lim  supjTy  _  -  t|  +  0  as  e  +  0  . 

T-k»  E1,ur 


But  (Oj,)  requires  that 

lim  sup)Tu£T  -  UKu^.)!  4  0  as  e  -*•  0  , 

f-Ko  9  T 

from  which  it  follows  simply  that  DKUj.)  +  t,  as  required.  □ 

Theorems  4.2  and  4.6  may  be  related  to  the  corresponding  results  for 
i.i.d.  sequences  in  the  following  way. 

Theorem  4.7.  Let  {u^,}  be  a  fcanily  of  constants  such  that  the  conditions 
of  Theorem  4.6  holdt  let  0  <  p  <  1*  and  let  h  be  chosen  as  in  (4.1) 
and  (4.8).  Then 

(4.15)  P{M(T)suT)  4  p  as  T  +  m 

if  and  only  if  there  is  a  sequence  (cn>  of  i.i.d.  random  variables  with 
common  d.f.  F  satisfying  1  -  F(u)  ~  hi{i(u)  as  u  -*■  <*>  and  such  that 
Mn  ■  ■«*(C1,?2...Cn)  eatisfies 
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(4.16)  P{VUnh}  p  * 

Proof.  If  there  is  an  i.i.d.  sequence  {cn}  with  common  d.f.  F  such  that 

(4.16)  holds  then  (as  noted  in  the  introduction)  we  have  1  -  F(unk)  ~  x/n, 
where  p  =  e‘T .  Since  1  -  F(u)  ~  hty(u)  we  have  iKu^)  ~  x/nh,  from  which 
(by  (4.8))  ijr(u-)  ~  T/T.  Hence  Theorem  4.2  gives  PtM^)^^}  •*  e'T  so  that 
(4.15)  holds. 

Conversely  if  (4.15)  holds  it  follows  from  Theorem  4.6  that 
TiKOj,)  -*■  T  and  hence  nto|/(unh)  x.  Let  (?n>  be  i.i.d.  random  variables 
with  the  same  d.f.  F,  say,  as  M(h),  so  that  by  (4.1) 

1  "  F(unh5  ~  ‘‘W'W  ~  T/n  * 

from  which  it  follows  that  Mn  ■  maxUj,^*  •  .Cn)  satisfies 

P(M  su  )  e"T  =  p,  as  required.  □ 

n  n 

These  results  show  how  the  function  may  be  used  in  the  classical 
criteria  for  domains  of  attraction  to  determine  the  asymptotic  distribution 
of  M(T) .  We  write  P(G)  for  the  domain  of  attraction  to  the  (extreme  value) 
d.f.  G,  i.e.  the  set  of  all  d.f.'s  F  such  that  Fn(x/an  +  bn)  -*■  G(x)  for 
some  sequences  {an>0),bn. 

Theorem  4.8.  Suppose  that  the  oonditione  of  Theorem  4.6  hold  for  all 
fmiliee  uJ  ■  x/aT  ♦  bT,  -«  <  x  <  «,  when  {aT>0),{bT)  are  given 
oonetants  and 

(4.17)  P{aT(M(T)-bT)  s  x)  G(x)  . 


Then 
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(4.18)  i|»(u)  ~  1  -  F(u)  as  u  -*■  «°  for  some  F  e  V(G )  . 

Conversely  if  (4.1)  holds  and  ij>(u)  satisfies  (4.18)  there  are  families  of 
constants  {a^,>0} , {b^,}  such  that  (4.17)  holds ,  provided  that  the  conditions 
of  Theorem  4.6  are  satisfied  for  each  -  x/a^  +  -<*>  <  x  <  ®. 

Proof.  If  (4.17)  holds,  together  with  the  conditions  stated.  Theorem  4.7 
shows  that 


H»nlA- bnh>  S  *>  -  GW 

where  M  is  the  maximum  of  n  i.i.d.  random  variables  with  a  common  d.f.  F_ , 
n  O’ 

say,  and  where  hij>(u)  ~  1  -  FQ(u)  as  u  an<*  fq  6  ’  We  may  choose 

a  d.f.  F  such  that  1  -  F(u)  *  ^(1-F0(u))  when  u  is  large  and  the  classical 
domain  of  attraction  criteria  show  that  F  e  0(G).  But  i|>(u)  ~  1  -  F(u)  as 
desired,  showing  (4.18). 

Conversely  if  (4.18)  holds  and  h  >  0  we  may  choose  FQ  e  D(G)  such  that 
hij>(u)  ~  1  -  FQ(u)  and  hence  define  an  i.i.d.  sequence  (Cn)  with  common 
d.f.  Fq,  An  =  max^,^.  •  *Cn) »  such  that 

P{a^(Mn-b^)  S  x)  -*■  G(x) 

for  some  constants  a^>0  ,  b^  .  Define  =  a^  ,  b^.*b^  for  nh  s  T  <  (n+l)h, 
n  =  0,1,2....  Then  (4.16)  holds  with  p  ■  G(x).  If  the  conditions  of 


Theorem  4.6  hold  for  each  vtj,  *  x/Sj  *  bT  then  (4.15)  holds,  which  yields 
(4.17).  □ 
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5.  Particular  classes  of  processes. 

In  this  section  we  first  show  how  the  conditions  required  for  the 
previous  theory  may  be  simplified  when  the  mean  number  y(u)  of  upcrossings 
of  each  level  u  by  £(t)  per  unit  time  is  finite,  and  then  briefly  indicate 
applications  to  stationary  normal  processes  (whether  or  not  y(u)  <  » ) . 
Throughout  Nu(I)  (Ny(t))  will  denote  the  number  of  upcrossings  of  the  level 
u  in  the  interval  I  (or  in  (0,t)  respectively)  . 

First  we  write  for  q  >  0 

(5.1)  Iq(u)  -  PU(0)<u<£(q)}/q  . 

Clearly  I  (u)  i  P{N  (q)*l}/q  S  EN  (q)/q  ■  V.  Further,  it  is  readily 
H  u  “ 

shown  (by  a  standard  dissection  of  the  unit  interval  into  subintervals  of 
length  q)  that 

(5.2)  y(u)  ■  lim  I  (u)  , 

q-*0  4 

which,  for  now,  we  assume  finite  for  each  u.  It  is  apparent  from  (5.2)  that 
y(u)  may,  at  least  in  principle,  be  readily  calculated  from  the  bivariate 
distributions  of  the  process.  It  may  also  happen  (as  for  many  normal 
processes)  that  I^(u)  ~  y(u)  as  u+«  when  q  depends  on  u,  q  »  q(u)  -►  0. 
For  greater  flexibility  we  shall  use  the  following  variant  of  such  a 
property.  Specifically  we  shall  assume,  when  needed,  that  for  each  a  >  0 
there  is  a  family  {q  (u)  -*■  0  as  u  -*■  »}  such  that  (with  q  ■  q  (u)  ,  y  ■y(u)) 

(5.3)  lim  inf  I  (u)/y  i  v 

u-*»  qa  8 

where  v(  +  1  as  a  +  0.  As  indicated  below,  for  many  normal  processes  we 
may  take  q#(u)  ■  a/u  and  snore  generally  as  aP{£(0)>u}/y(u)  . 


Jtmi 
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We  shall  assume  as  needed  that 

(5.4)  P{£(0)>u}  ■  oy(u)  as  u  , 


which  holds  under  general  conditions.  For  example,  it  is  readily  verified 
if  for  some  q  =  q(u)  -*■  0  as  u  -*•  «, 


(5.5) 


lim  sup 

U-H» 


<  1 


since  (5.5)  implies  that  lim  inf  ql  (u)/P{£(0)>u)  >  0,  from  which  it 
follows  that  P{£(0)>u}/I^(u)  0,  and  hence  (5.4)  holds  since 

Iq(u)  £  y(u) . 

We  may  now  recast  the  conditions  (3.6)  and  (3.7)  in  terms  of  the 
function  y(u). 


Lemma  5.1. 

(i)  Suppose  y(u)  <  00  for  each  u  and  that  (5.4)  (or  the  sufficient 
condition  (5.5))  holds.  Then  (3.6)  holds  with  <|>(u)  =  y(u)  . 
(ii)  If  (5.3)  holds  (for  some  family  {q  (u)})  then  (3.7)  holds  with 

ft 

tfi(u)  *  P(u)  . 

Proof.  Since  clearly 

P{M(h)>u)  s  P{Nu(h)il}  +  P(5(0)>u)  s  yh  +  P{£(0)>u}  , 


(3.6)  follows  at  once  from  (5.4),  which  proves  (i) . 

How  if  (5.3)  holds,  then  with  q*qg(u)  ,  y»y(u), 

P(C(0)<u,  C(q)<u,M(q)>u)  -  P(C(0)<u  ,  M(q)>u)  -  PU(0)<u<  C(q)> 

S  P{Nu(q)il)  -  qlq(u) 

S  yq  -  yqva(l*o(l)) 
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so  that 


lim  sup  P(C(0)<u  ,  £(q)<u  ,  M(q)>u}/(qy)  s  1  -  v  , 

xr+oo  8 

which  tends  to  zero  as  a  0,  giving  (3.7).  □ 

In  view  of  this  lemma,  Gnedenko's  Theorem  now  applies  to  processes  of 
this  kind  using  the  more  readily  verifiable  conditions  (5.3)  and  (5.4),  as 
follows. 

Theorem  S.2.  Theorem  3.5  holds  for  a  stationary  prooess  £(t)  with 

tp(u)  «  jj(u)  <  “  for  each  u  if  the  conditions  (3.6)  and  (3.7)  are  replaced 

by  (5.4)  and  (5.3).  □ 

Finally,  the  condition  D^(Uj.)  may,  in  certain  circumstances,  be 
replaced  by  a  sufficient  condition  involving  the  second  moment  of  Nu(l)  when 
this  is  finite.  This  condition  is  not  necessarily  simpler  to  verify,  but 
the  second  moment  involved  may  usually  be  obtained  in  terms  of  (integrals 
containing)  the  joint  densities  of  the  process  and  its  derivative  at  two 
general  points  tj , tj  . 

2 

Lemma  5.3.  Suppose  that  for  the  stationary  prooess  £(t),  EN‘(1)  <  »,  and 
for  a  given  family  (u)  ■  {u^,} , 

(5.6)  i  lim  sup  EN11(eT)(N„(eT)-l)  «  e  +  0. 

E  T-*«  u  u 


Then  £(t)  satisfies  D*^),  with  ♦(u)  ■  p(u)  »  ENu(l)  . 
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Proof.  Clearly,  writing  y  ■  wCu^.) , 

°  S  T(l'-U£T,U)  ■  £  E(Nu<CT)-N€T.U'£T»  • 

Now  if  Nu(eT)  >  1,  Nu(eT)  -  N^^CcT)  s  Nu(£T) (NJeTj-l) .  Also 
Nu(eT)  -  N£T  u(eT)  is  zero  if  Nu(eT)  ■  0  and  is  zero  or  1  if  Nu(eT)  *  1, 
the  latter  case  requiring  that  Nu(-eT,0)  Js  1  also.  Hence  we  have 

0  s  E{Nu(eT)  -  NeT>u(eT)}  s  ENu(eT) (NU(ET)-1)  +  P{N(2eT)>1> 
s  ENu(£T)(Nu(eT)-1)  ♦  ENu(2eT)(Nu(2eT)-1) 

so  that  Dc(tUj.)  follows  by  applying  (5.6)  twice  (once  with  2e  replacing  c) .  □ 

For  stationary  normal  processes,  finiteness  of  EN^fl)  requires  a  little 
more  than  existence  of  the  second  spectral  moment  used  to  ensure  finiteness 
of  jj  (cf.  [3]).  We  turn  now  to  the  consideration  of  stationary  normal 
processes,  but  will  not  restrict  attention  to  those  for  which  even 
y  •  EN 4(1)  is  finite.  Specifically  we  assume  that  £(t)  is  a  (zero  mean) 
stationary  normal  process  with  covariance  function 

(5.7)  r(r)  »  1  -  C|t|°  o|t|q  u  t  +  0 

for  some  a,  0  <  a  s  2.  (The  case  o  ■  2  gives  y  <  • .)  There  is  a 
considerable  literature  dealing  with  extremal  properties  of  such  processes, 
and  of  slightly  more  general  cases  (which  could  be  included  here)  in  which 
the  term  |r|a  is  multiplied  by  a  slowly  varying  function  as  t  ♦  0 
(cf.  [2],  [16]).  Of  course  a  number  of  the  sa me  arguments  (which  in  some 
cases  are  rather  intricate)  used  in  these  papers  are  required  to  verify  our 
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general  conditions  here.  We  will  not  attempt  to  reproduce  these  arguments 
but  rather  to  simply  indicate  the  basic  considerations  used  and  where  they 
may  be  found.  However  it  will  be  convenient  to  summarize  these  results  as  a 
theorem  even  though  formal  proofs  are  not  given. 


Theorem  5.4.  Let  £(t)  be  a  zero  mean  stationary  normal  process  with 
covariance  function  r(t)  satiefying  (5.7).  Then 

(i)  (3.6),  and  in  fact  (4.1),  hold  with  t|»(u)  •  C1/ot  Hau2/a  <|>(u)/u, 
in  which  f  is  the  standard  normal  density  C  is  as  in  (5.7),  and 
Ha  is  a  constant  depending  only  on  a . 

(ii)  (3.7)  holds  with  q&(u)  ■  su~2^a  . 

(Hi)  Dc(Oj.)  holds  with  respect  to  a  family  (q)  if  T4»(uuj.)  is  bounded 
and 


(5.8) 


I 

q 


l 

XT*  kq  S  T 


|r(kq)|e 


-u£/(Hr(kq)|) 


0  as  T  ♦  ® 


for  each  X  >  0 .  This  holds ,  in  particular ,  if  TiK^)  is 
bounded  (with  iji  defined  as  in  (i))  and  r(t)  log  t  0  as  t  -*•  ®. 
(iv)  If  r(t)  log  t  *  0  and  HK*^)  ■*>  T  >  0,  then  D'O^)  holds  and 
P(M(T)s  u?}  -►  e'T. 

(v)  If  r(t)  log  t  ♦  0,  M(T)  has  the  limiting  distribution  given  by 


P(aT(M(T)-bT)  s  x)  e"e 


fc,.  -  (2  log  T)% 

bT  -  (2  log  T)**  ♦  (2  log  T)”**  (( j-y)10g  log  T 
♦  10g(2°  1-1  W~**  C1/a  Ha)>  . 


where 


Indications  and  sources  of  proof. 

(i)  A  derivation  of  (4.1)  (froa  which  (3.6)  follows)  appears  in 

several  developments  of  the  normal  theory  (e.g.  Theorem  2.1  of  [16]).  In 

the  case  a  ■  2,  (3.6)  is  incidentally  simply  obtained  froa  "Rice's  formula" 

2 

V  “  [(-r"(o))%/2ir]e‘u  /2. 

(ii)  This  may  be  shown,  for  example,  along  the  lines  of  Lemma  2.4  of 
[16] ,  although  a  more  direct  derivation  is  obtainable  from  the  normal  theory 
given  by  Lindgren  and  Rootzen  in  [13]. 

(iii)  The  proof  of  this  involves  a  standard  calculation  using  "Slepian's 
Lemma"  (cf.  Lemma  3.5  of  [15]),  froa  which  it  follows  that  for  two  sets  of 
standard  normal  random  variables  ^  with  covariance  matrices 

l*ij] > [vijl  »  lXijMyijl 


n  n  ,  ,  -u2/(l+|X, . 

|P{  n  (?.su)}  -  P{  n  (n,su)}|  s  K  l  |X..  -  v..|(l-X2  )-J*  e  ij 

j«l  J  j«l  J 


I) 


In  this  application  (using  the  notation  of  (3.4)),  the  ^  are  identified 
with  the  r.v.'s  C(Sj)...  5(*p)  •  ^(tj)...  C(tpf)  and  the  with  p+p' 
standard  normal  r.v.'s  having  the  same  correlations  except  that 
cov(C(si) ,  C(tj))  is  replaced  by  zero  for  isisp  ,  lsjsp'. 

The  fact  that  boundedness  of  T9(UjJ  together  with  r(t)  log  t  *  0 
implies  (5.8)  follows  by  standard  calculations  (cf.  [1]  or  Lemma  3.1  of  [13]) 

(iv)  If  r(t)  log  t  ♦  0  and  Tty(Uj.)  x  >  0  then  D^(ttj.)  may  be 
obtained  from  arguments  leading  to  Theorem  3.1  of  [15],  though  it  seems 
likely  that  a  shorter  route  via  our  Lemma  3.3  may  be  possible.  It  then 
follows  from  Theorem  4.2  that  P{M(T)f  uT>  -*>  e"T.  Of  course  any  proof  (of 
which  there  are  several)  that  r(t)  log  t  ♦  0  and  Tf^)  t  implies  that 
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PCmOOsu^J  +  e_T  must  also  imply  that  D^fu^.)  holds  by  virtue  of  our 
]  n— r  4.5.  That  is  D'(Uj.)  may  be  regarded  as  a  necessary  condition  for 
(4.3)  to  imply  (4.4). 

(v)  This  follows  at  once  from  the  (relatively)  straightforward 
verification  of  the  fact  that  Tty(u.j.)  -►  T  *  e  x  when  u^,  «  x/fcj.  ♦  b.j. ,  using 
the  above  results.  D 

6.  Poisson  and  related  properties. 

In  this  section  we  shall  just  briefly  indicate  the  Poisson  properties  j 

associated  with  high  level  upcrossings.  We  confine  the  discussion  to  the  j 

case  where  the  number  Nu(I)  of  upcrossings  in  a  bounded  interval  I  has  a  I 

finite  mean,  writing  again  y  ■  y(u)  *  ENu(1).  Cases  where  this  is  not  so  j 

J 

are  similarly  dealt  with  in  terms  of  e-upcrossings. 

Our  objective  is  to  show,  under  D£  and  conditions,  that  the  point  j 

process  of  upcrossings  of  a  high  level  takes  on  a  Poisson  character — as  is 
well-known  in  the  case  when  the  stationary  process  £(t)  is  normal.  Since 
the  upcrossings  of  increasingly  high  levels  will  tend  to  become  rare,  a 
normalization  is  required.  To  that  end  we  consider  a  time  period  T  and  a 
level  uT,  both  increasing  in  such  a  way  that  Ty  -*-t  ,  (y  *  y(u^.)) ,  and 
define  a  normalized  point  process  of  upcrossings  by 

-N(tT)) 

HP 

for  each  interval  (or  more  general  Borel  set)  I,  so  that,  in  particular. 


N*(I)  -  N^(TI)  ,  (N£(t) 


(6.1) 


EN*(1)  ■  EHj.fT)  ■  VT  ♦  T  . 


This  shows  that  the  "intensity"  (i.e.  nean  nuaber  of  events  per  unit 
tine)  of  the  (normalized)  upcrossing  point  process  converges  to  t.  Our  task 
is  to  show  that  the  upcrossing  point  process  actually  converges  (weakly)  to 
a  Poisson  process  with  nean  T. 

The  derivation  of  this  result  is  based  on  the  following  two  extensions 
of  Theorem  4.2,  which  are  proved  by  similar  arguments  to  those  used  in 
obtaining  Theorem  4.2. 

Theorem  6.1.  Under  the  conditions  of  Theorem  4.2 ,  if  0  <  1  and  pT  t, 
then 

(6.2)  •  P{M(0T)iuT}  «*•  e‘6T  fli  T  ♦  •  .  0 

Theorem  6.2.  If  Ij, I2- • . 1^  are  disjoint  subintervals  of  [0,1]  and 

IJ  »  TIj  ■  (t :  t/T  e  1}  then  under  the  conditions  of  Theorem  4.2t  if  pT  t, 

k  k 

(6.3)  P{  n  MapSu,.)  -  P{M(Ij)SuT) -*•  0  , 
so  that  by  Theorem  8.1 

k  -ire. 

(6.4)  P{  n  (W(IJ)s uT)}  -*>  e  , 

j-1  3  T 

ahere  8^  is  the  length  of  ly  l  s  )  s  k .  0 

It  is  now  a  relatively  straightforward  matter  to  show  that  the  point 
processes  N£  converge  (in  the  full  sense  of  weak  convergence)  to  a  Poisson 
process  N  with  intensity  t. 
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Theorem  6.3.  Under  the  conditions  of  Theorem  4.2 ,  if  Ty  ■*  x  where 
y  ■  y(u.j.),  then  the  family  N*  of  (normalised)  point  processes  of  uperossings 
of  Uj,  on  the  unit  interval  converges  in  distribution  to  a  Poisson  process  N 
with  intensity  x  on  the  unit  interval  as  T  «. 

Proof.  By  Theorem  4.7  of  [8]  it  is  sufficient  to  prove  that 

(i)  EN*{ (a,b] }  -*>  EN{(a,b]}  -  t(b-a)  as  T  «  for  all  a,b, 

0  s  a  s  b  s  1. 

n 

(ii)  P{N*(B)*0}  -*•  P(N(B)«0}  as  T  +  *  for  all  sets  B  of  the  fora  uB. 

1  1  1 
where  n  is  any  integer  and  B^  are  disjoint  intervals 

(a.,b4]  c  (0,1]  . 

Now  (i)  follows  trivially  since 

En*{ (a,b]}  «  yT(b-a)  T(b-a)  . 

To  obtain  (ii)  we  note  that 

0  s  P{N*(B)-0}  -  P{M(TB)sut> 

-  P{Nu(TB)-0  ,  M(TB)>uT) 

*  l  P  (5(Ta.  )>u_} 
i-1  1  1 

n 

since  if  the  maximum  in  TB  ■  u  (Ta.,Tb,]  exceeds  u_,  but  there  are  no 

i-1  1  1  1 

uperossings  of  Uj.  in  these  intervals,  then  €  must  exceed  u  at  the  initial 
point  of  at  least  one  such  interval.  But  the  last  expression  is  just 
nP{C(0)>Uj.}  o  as  T  ■*  •.  Hence 
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P(H*(B)-0}  -  P{M(TB)£ Uy)  -*•  0  . 


-  a.f 

But  POfCrB)*!^}  •  P{  n  (M(TBj)< ttj.}  ♦  e'**  by  Theorem  6.2  so  that 

. . . - =n:£(b.-  a.)  -  _ 

(ii)  follows  since  P{N(B)-0}  »  e  1  □ 


Corollary.  If  B^  are  disjoint  (Borel)  subsets  of  the  unit  interval  and 
if  the  boundary  of  each  IL  hae  aero  Lebeegue  measure  then 


PttfCB^-rj  , 


*1 

n  -Tm(B.)  [T»(B.)J 
lsisn)  -►  n  e  1  - i - 

i-1  V 


where  m(B..)  denotes  the  Lebeegue  measure  of  B^ . 


Proof.  This  is  an  immediate  consequence  of  the  full  weak  convergence  proved 
(cf.  Lewi  4.4  of  [8]).  □ 

The  above  results  concern  convergence  of  the  point  processes  of 
upcrossings  of  Uj.  in  the  unit  interval  to  a  Poisson  process  in  the  unit 
interval.  A  slight  modification  (requiring  D  and  D'  to  hold  for  all 

V  C 

families  uflT  in  place  of  u^,  for  all  8  >  0)  enables  a  corresponding  result 
to  be  shown  for  the  upcrossings  on  the  whole  positive  real  line,  but  we  do 
not  pursue  this  here.  Instead  we  show  how  Theorem  6.3  yields  the  asymptotic 
distribution  of  the  rth  largest  local  maximum  in  (0,T). 

Suppose,  then,  that  £(t)  has  a  continuous  derivative  a.s.  and  define 
N^(T)  to  be  the  number  of  local  maxima  in  the  interval  (0,T)  for  which  the 
process  value  exceeds  u,  i.e.  the  number  of  downcrossing  points  t  of  zero  by 
C  in  (0,T)  such  that  C(t)  >  u.  Clearly  N^(T)  a  NU(T)  -  1  since  at  least 
one  local  maximum  occurs  between  two  upcrossings.  It  is  also  reasonable  to 
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expect  that  if  the  staple  function  behavior  is  not  too  irregular  that  there 
will  tend  to  be  just  one  local  naxiaua  between  most  successive  upcrossings 
of  u  when  u  is  large,  so  that  N^(T)  and  NU(T)  will  tend  to  be  approximately 
equal.  The  following  result  makes  this  precise. 

Theorem  6.4.  With  the  above  notation  let  {u^}  be  oomtante  euoh  that 
TyC-TyfUj))  t  >  0.  Suppose  that  En^(1)  is  finite  for  each  u  and  that 
ENy(l)  ~  y(u)  ae  u  «.  Then,  writing  Oj.  ■  u,  E|N^(T)  -  Nu(T)  |  0. 

If  also  the  conditions  of  Theorem  6.3  hold  (so  that 
P{Nu(T)-r}  e'V/r!)  it  follows  that  P(N^(T)-r}  e'V/r!  . 

Proof.  As  noted  above,  N^(T)  2  NU(T)  “  ^  it  is  clear,  moreover,  that 
if  N^(T)  -  Nu(T)  -  1,  then  C(T)  >  u.  Hence 

e|N^(T)  -  Nu(T)|  -  E{N^(T)  -  Nu(T)}  ♦  2P{N^(T)  -  NJT)  -  1} 

S  TEN^(l)  -  yT  ♦  2P{£(T)>u}  , 

which  tends  to  zero  as  T  ■*  •  since  P{?(T)>u^}  ■  P{?(0)>uT}  -►  0  and 

TEN'  (1)  -  yT  ■  yT[(l+o(l))  -  1]  0,  so  that  the  first  part  of  the  theorem 

UT 

follows.  The  second  part  now  follows  immediately  since  the  integer-valued 
r.v.  N^(T)  -  NU(T)  tends  to  zero  in  probability,  giving  P(N^(T)  «  NU(T)}  -*•  0 
and  hence  P{N^(T)«r}  -  P{Nu(T)»r}  -*■  0  for  each  r.  □ 

Now  write  M^(T)  for  the  rtJ>  largest  local  maximum  in  the  interval 
(0,T).  Since  the  events  (T)su},{N^(T)<  r}  are  identical  we  obtain  the 
following  corollary: 


Corollary  1.  Undo*  the  oonditiona  of  the  theorem 

P{M(r)(T)su_}  -►  e"T  1  xVsl  .  □ 

1  s-1 

As  a  further  corollary  we  obtain  the  Uniting  distribution  of  (T) 
in  terms  of  that  for  M(T) . 

Corollary  2.  Suppose  that  P{aT(M(T)-bT)  £  x}  ■*>  G(x)  and  that  the 

oonditiona  of  Theorem  4.6  hold  with  Uj.  ■  X/*T  +  bT  for  each  real  x  (and  ip  »  y). 

Suppoee  also  that  En^(1)  ~  ENu(l)  as  u  -*•  ®.  Then 

(6.5)  P{aT(M(r)(T)-bT)  S  x}-*-  G(x)  V  I-lo8  G(x))s/sl  , 
where  G(x)  >  0  (and  zero  if  G(x)  ■  0). 

Proof.  This  follows  from  Corollary  1  by  writing  G(x)  ■  e"T  since 
Theoren  4.6  implies  that  Ty  -*■  t  .  □ 

Note  that  for  a  stationary  normal  process  with  finite  second  and  fourth 
spectral  moments  X2,  it  may  be  shown  (Section  11.6  of  [3])  that 

OT(l)  «  y*(uX2/d)  ♦  (A4/X2),i  [1  -  *{u(A4/A),s} 

2 

where  A  »  X4  -  X2  and  4  is  the  standard  normal  d.f . ,  so  that  clearly 
EN^(l)  ~  y  as  u  • . 

The  relation  (6.S)  gives  the  asymptotic  distribution  of  the  rth  largest 
local  maximum  M^(T)  as  a  corollary  of  the  Poisson  result.  Theorem  6.4. 

This  Poisson  result  may  itself  be  generalized  to  apply  to  joint  convergence 


38 


of  upcrossings  of  several  levels  to  a  point  process  in  the  plane  composed  of 
successive  "thinnings"  of  a  Poisson  process.  From  a  result  of  this  kind  it 
is  possible  to  obtain  the  joint  asymptotic  distribution  of  any  number  of  the 
M^r) (T) ,  and  also  of  their  time  locations. 
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